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Abstract 


voltage  and  currents  are  changing  in  time,  the  impedance 


A  generalized  magnetically  insulated  transmission  line 
(MITL)  flow  model  has  been  developed  to  treat  dynamic 
MITL  problems  [1].  By  including  electron  pressure  in  the 
model  and  allowing  non-zero  values  of  the  electric  field  at 
the  cathode,  this  MITL  model  can  treat  both  emission  and 
re-trapping  of  flow  electrons.  Most  previous  MITL  flow 
models  only  describe  equilibrium  flow  conditions  without 
emission  or  re-trapping  and  cannot  adequately  treat 
dynamic  situations.  Such  dynamic  situations  are  common 
and  include  impedance  transitions  along  the  line,  variable 
impedance  transmission  lines,  coupling  to  loads,  etc.,  all 
of  which  can  cause  electron  emission  from  the  cathode 
and/or  electron  re-trapping  onto  the  cathode.  The  model 
is  being  benchmarked  against  particle-in-cell  (PIC) 
simulations  using  the  LSP  code  [2].  Of  particular  interest 
for  this  benchmarking  effort  is  the  treatment  of  re¬ 
trapping  waves  that  occur  when  the  MITL  is  terminated 
by  an  under-matched  load.  Ultimately,  the  model  will  be 
incorporated  into  a  transmission  line  code  such  a 
BERTHA  [3]  so  that  MITL  problems  can  be  studied  more 
quickly  and  efficiently  than  with  PIC  codes. 


I.  MITL  MODEL 

Many  modem  pulsed  power  generators  use  magnetically 
insulated  transmission  lines  (MITL)  to  couple  power 
between  the  driver  and  the  load.  In  an  MITL  the  electric 
field  stress  on  the  cathode  exceeds  the  vacuum  explosive- 
emission  threshold  and  electron  emission  occurs.  For 
sufficiently  high  current,  emitted  electrons  are 
magnetically  insulated  from  crossing  the  anode-cathode 
gap  and  flow  axially  downstream  in  the  direction  of 
power  flow  as  illustrated  in  Fig.  1.  The  return  current 
from  the  total  anode  current  L  is  divided  between  current 
Ic  flowing  in  the  cathode  and  current  If  flowing  in  the 
vacuum  electron  layer,  i.e..  If  =  L  -  L-  As  a  result  of  the 
electron  flow  in  vacuum  between  the  electrodes  the 
impedance  of  the  MITL  is  altered  and,  thus,  the  power 
coupling  between  the  machine  and  the  load  changes.  For 
equilibrium  flow  it  has  been  shown  that  the  effective 
impedance  of  the  MITL  is  best  described  by  the  flow 
impedance  Zf  [4,5].  In  a  dynamic  system  where  the 
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Figure  1.  Schematic  of  MITL  flow  in  negative  polarity. 


also  varies  in  time  along  the  line  [1,6]. 

A  new  generalized  model  for  MITL  flow  was 
developed  in  [1]  for  incorporation  into  the  transmission 
line  code  (TLC)  Bertha  [3]  to  treat  dynamic  MITL 
problems  so  that  a  more  computationally  intensive  PIC 
code  treatment  is  not  required.  The  model  describes  both 
self-limited  flow  as  the  pulse  initially  propagates  down 
the  MITL  toward  the  load  and  the  subsequent  electron 
power  flow  along  the  MITL  after  the  pulse  encounters  the 
load.  To  accomplish  this,  the  model  must  treat  electron 
emission  at  the  pulse  front  as  illustrated  in  Fig.  2  and  at 
impedance  transitions  along  the  MITL  where  required. 
For  low  impedance  loads,  this  description  also  includes 
electron  re -trapping  [7]  as  the  flow  is  modified  by  the 
wave  reflection  off  the  load  and  the  percentage  of  the 
return  current  in  vacuum  electron  flow  decreases  as 
illustrated  in  Fig.  3. 

The  model  provides  four  generalized  MITL  flow 
equations  that  must  be  solved  simultaneously  and  are 
presented  as  Eqs.  (22)  -  (28)  in  [1].  For  given  values  of 
V,  L,  and  A,  these  equations  can  be  solved  for  Qa,  Q^,  L, 


Figure  2.  Schematic  of  schematic  of  self-limited  flow 
where  L  >  0  implies  electron  emission  from  the 
cathode,  L  <  0  implies  electron  re -trapping  onto  the 
cathode,  and  J»  >0  imnlies  electron  loss  to  the  anode. 
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Figure  3.  Schematic  showing  a  re-trapping  wave 
generated  by  the  interaction  of  the  MITL  flow  with  a 
load  impedance  that  is  less  than  the  self-limited 
impedance  of  the  line. 

and  Zf,  where  Qa  and  Qc  are  the  charge  per  unit  length  at 
the  anode  and  cathode  respectively.  The  MITL  is 
modeled  by  a  series  of  one-time-step-long  transmission 
line  elements  characterized  by  their  capacitive  impedance 
Zcap  and  a  loss  resistor  Rioss  as  illustrated  in  Fig.  4.  For  a 


time  step  of  x,  the  element  has  length  f  =  ex  where  c  is  the 
speed  of  light.  Relating  fQa  =  CV  =  xV/Zcap,  the 
capacitive  impedance  is  defined  in  [1]  in  terms  of  Zf  as 

y  _  ^O^Qc  +^f(^Qa  ~^Qc) 


cap 
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Figure  4.  Representation  of  an  MITL  TLC  element 
showing  forward-  and  backward-going  voltage  waves  at 
each  end  of  the  element. 


ZoJc=- 


5(Z„cQj  5(ZjJ 


(3) 


Note  that  in  the  basic  equilibrium  MITL  model  Qc  =  0  so 
that  in  this  special  case  Zcap  =  Zf.  The  loss  resistor 
represents  any  electron  current  loss  to  the  anode  in  an 
element  and  is  defined  below. 

The  values  for  V  and  la  are  obtained  as  a  function  of 
time  by  advancing  the  TLC  in  time  and  the  values  of  A 
are  obtained  by  appropriately  solving 
d(cA)  _  ^ 

dt  dz 

(2) 

at  each  element  for  each  time  step.  Between  time  steps 
the  Zcap  of  each  element  is  adjusted  to  reflect  the  new  flow 
properties  of  the  element.  In  order  to  preserve  charge  and 
flux  when  the  element  impedance  is  changed  the  waves 
on  the  TLC  element  must  also  be  adjusted. 

The  linear  electron  current  density  Jc  associated  with 
emission  from  the  cathode  or  re-trapping  onto  the  cathode 
is  determined  by 


cdt  dz 

Here,  positive  Jc  implies  electron  flow  in  the  direction 
from  the  cathode  to  the  anode  in  negative  polarity  so  that 
Jc  >  0  implies  emission  and  Jc<0  implies  re-trapping. 
Similarly,  electron  current  loss  to  the  anode  boss  =  f  Ja  is 
described  by 

7  ,  dZpCQ.)  3(z,l.) 

°  '  ca  dz  ' 

Once  boss  is  found,  Rioss  =  V/boss  is  applied  at  the  next  time 
step. 

A  sample  set  of  solutions  for  V  =  6  MV  and  cA  =  1 
MV  as  a  function  of  la  (shown  as  the  vertical  axis)  is 
given  in  Figs.  5-9.  In  these  figures  valid  solutions  are 
restricted  to  the  clear  areas  and  the  dashed  line  indicates 
the  location  of  the  Qc  =  0  basic  equilibrium  solutions  for 
V  =  6  MV;  however,  the  value  for  cA  varies  along  this 
line  and  is  not  fixed.  The  cA  =  7  MV  case  for  the  basic 
solution  is  indicated  by  the  open  circle  where  the 
generalized  solution  and  the  basic  solution  curves  cross 
(see  [1]  for  more  details).  The  open  squares  signify  the 
two  saturated  flow  solutions  of  the  generalized  flow 
model  as  derived  in  [1]. 


Figure  5.  Plot  of  Qa  vs.  Q  solutions. 
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Figure  6.  Plot  of  Qc  vs.  la  solutions. 
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Figure  7.  Plot  of  Ic  vs.  la  solutions. 


Figure  9.  Plot  of  Zcap  vs.  la  solutions. 


Before  the  transmission  line  begins  to  emit,  vacuum 
solutions  with  no  electron  flow  apply  with  ZocQa  =  ZqcQc 
=  V,  Zola  =  ZqIc  =  cA,  and  Zf  =  Zq.  As  an  example,  the 
vacuum  solution  is  indicated  in  Fig.  7.  Once  the  electric 
field  exceeds  the  explosive  emission  threshold  and  the 
cathode  is  turned  on  to  emission,  MITL  flow  solutions  are 
required.  For  the  case  shown  in  Fig.  7  with  V  =  6  MV, 
cA  =  7  MV  and  a  given  value  of  la,  there  are  three 


possible  outcomes:  one  solution  is  found  (upper  left  of  the 
curve),  two  solutions  are  found  (lower  right  of  the  curve), 
or  no  solution  is  found  (between  vacuum  solution  and 
curve).  The  special  solution  indicated  by  the  line 
connecting  the  vacuum  solution  and  the  saturated  flow 
solution  is  used  here  to  model  the  transition  between  the 
vacuum  line  and  the  MITL  flow  case. 

II.  BENCHMARKING  MODEL  RESULTS 
WITH  PIC  SIMULATIONS 

Although  more  benchmarking  is  needed,  an  initial  test 
of  the  generalized  MITL  flow  model  is  provided  by 
compared  theoretical  predictions  from  the  model  with  a 
simple  nearly  steady  state  simulation.  The  simulation  was 
run  using  the  LSP  code  [2]  for  an  MITL  with  a  vacuum 
impedance  of  30  Emission  is  turned  off  in  the  diode 
load  region  so  that  it  presents  an  open  circuit  to  the  MITL 
and  the  flow  remains  self-limited  as  illustrated  in  Fig.  10. 
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Figure  10.  LSP  simulation  of  self-limited  flow. 

Properties  of  the  MITL  flow  at  the  location  indicated  as 
pos#3  in  Fig.  10  are  summarized  in  Table  1  along  with 
predictions  from  the  generalized  MITL  flow  model 
described  herein.  The  input  for  the  model  are  V  =  5.92 
MV,  cA  =  6.92  MV  and  Zola  =8.1  MV  and  are  presented 
in  red  in  the  table.  The  simulation  parameters  V  and  cA 
are  close  to  those  used  to  generate  Figs.  5-9  so  that  one 
can  see  roughly  where  this  specific  solution  lies  on  the 
solution  curves.  The  LSP  results  are  close  to  the  Qc  =  0 
solution  while  the  model  predicts  the  solution  to  lie 
slightly  to  the  left  of  it.  This  small  discrepancy  probably 
results  from  the  difference  between  the  density  profile  in 


Table  1.  MITL  flow  properties. 


LSP 

Model 

V(MV) 

5.92 

5.92 

cA(MV) 

6.92 

6.92 

Z„Ia(MV) 

8.10 

8.10 

ZoIc(MV) 

2.97 

3.42 

ZocQa(MV) 

7.50 

7.47 

ZocQc(MV) 

0.18 

-1.44 

Znan(O) 

23.7* 

23.7 

‘calculated  from  V/cQ^ 
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the  electron  layer  assumed  in  the  model  compared  with 
that  obtained  in  the  simulation.  This  discrepancy  mostly 
affects  Ic  and  Qc.  However,  the  result  that  affects  the  TLC 
description  of  the  power  flow  most  is  Z^ap  and  that 
prediction  is  very  accurate  as  is  Qa  which  is  related  to  Zcap 
through  Zcap  =  V/cQa. 

III.  INCORPORATING  THE  MITL 
MODEL  INTO  BERTHA 

Because  electron  space-charge  and  current  are 
distributed  in  the  line,  there  is  not,  in  general,  a  well- 
defined  single  wave  impedance  that  describes  an  MITL. 
Previously,  the  electric  flow  impedance  Zf  and  magnetic 
flow  impedances  Z^  were  defined  by  the  distances  of  the 
centroid  of  the  charge  and  the  centroid  of  the  current 
from  the  anode,  respectively  [4,6].  Recently,  the 
capacitive  impedance  Zcap  and  the  inductive  impedance 
Zind  were  introduced  in  the  generalized  MITL  flow  model 
and  related  to  Zf  and  Z^,  to  describe  the  electrical 
properties  of  the  MITL  [1].  In  fact,  Zcap  ^  Zf  for  the  basic 
model  where  ZqcQc  =  0.  The  difference  between  Zcap  and 
Zind  is  usually  small  so  that  Zcap  provides  a  reasonable 
approximation  of  the  MITL  flow  impedance  under  most 
circumstances  of  interest.  For  purposes  of  incorporating 
the  generalized  MITL  flow  model  into  a  TLC  where  a 
single  impedance  is  needed  to  characterize  a  transmission 
line  element,  Zcap  as  described  herein  is  used  to  describe 
the  MITL  impedance.  Alternately,  another  value  such  as 
the  average  value  (Zcap  +  could  be  used  to  describe 

the  MITL  impedance.  The  reader  is  referred  to  [1]  for  a 
description  of  Zjnd.  For  incorporating  the  model  into  a 
circuit  code  were  the  MITL  is  modeled  by  a  sequence  of 
series  inductors  and  parallel  capacitors,  both  Zcap  and  Zind 
could  be  used  to  describe  the  capacitance  and  inductance 
of  these  circuit  elements  respectively  to  provide  better 
simulation  fidelity. 

Because  it  is  assumed  that  electrons  react 
instantaneously  to  the  time-dependent  fields,  the 
generalized  MITL  flow  model  is  a  quasi-equilibrium 
model.  However,  when  combined  with  time-dependent 
circuit  equations  for  evolving  V,  L,  and  A  in  time,  this 
quasi-equilibrium  model  can  be  used  to  build  a  dynamic 
model  for  MITL  flow,  for  example  in  a  transmission  line 
code.  Work  is  now  underway  to  develop  a  robust  and  fast 
numerical  technique  for  solving  the  new  MITL  flow 
equations.  Ultimately  this  will  allow  efficient  and 
accurate  modeling  of  MITL  flow  in  a  fast  transmission 
line  code  to  replace  the  more  computationally  intensive 
particle-in-cell  code  treatment. 

IV.  SUMMARY 

A  generalized  magnetically  insulated  transmission  line 
(MITL)  flow  model  has  been  developed  to  treat  dynamic 
MITL  problems.  By  including  electron  pressure  in  the 
model  and  allowing  non-zero  values  of  the  electric  field  at 


the  cathode,  this  MITL  model  can  treat  both  emission  and 
re-trapping  of  flow  electrons.  Most  previous  MITL  flow 
models  only  describe  equilibrium  flow  conditions  without 
emission  or  re-trapping  and  cannot  adequately  treat 
dynamic  situations.  Such  dynamic  situations  are  common 
and  include  impedance  transitions  along  the  line,  variable 
impedance  transmission  lines,  coupling  to  loads,  etc.,  all 
of  which  can  cause  electron  emission  from  the  cathode 
and/or  electron  re-trapping  onto  the  cathode.  The  model 
is  now  being  benchmarked  against  particle-in-cell  (PIC) 
simulations.  The  initial  benchmarking  results  presented 
here  show  that  the  model  produces  a  very  accurate 
description  of  Zcap  which  is  essential  for  accurately 
advancing  a  TLC  in  time.  Of  particular  interest  for  this 
benchmarking  effort  is  the  treatment  of  re -trapping  waves 
that  occur  when  the  MITL  is  terminated  by  an  under¬ 
matched  load.  Future  work  will  focus  on  this  re-trapping 
phenomenon.  Ultimately,  the  model  will  be  incorporated 
into  a  transmission  line  code  such  a  BERTHA  so  that 
MITL  problems  can  be  studied  more  quickly  and 
efficiently  than  with  PIC  codes. 
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